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ABSTRACT 

In  the  present  paper,  the  equation  of  Magneto  hydrodynamic  boundary  layer  model  for  Newtonian  power  law 
jluid  has  been  analyied  in  the  presence  of  transverse  magnetic  field  for  two  cases,  (i)  when  plate  is  stationary  (ii)  when 
the  fluid  and  plate  moves  in  same  direction  and  with  same  velocity.  Here  we  solve  the  governing  nonlinear  differential 
equation  with  their  associated  boundary  conditions  using  Spline  functions  due  to  Blue.  The  beauty  of  this  method  is  we 
can  solve  nonlinear  problem  directly,  without  converting  to  the  linear  form.  Method  description  is  explained  and 
graphical  results  describing  the  displacement  and  velocity  of  the  flow  are  presented  in  this  work. 

KEYWORDS:  Power-Law  Fluids,  Magnetic  Field,  Nonlinear  Differential  Equation,  Collocation  Method  & Spline 
Functions 

Nomenclature: 

MHD  - Magneto  hydrodynamics 

- Yelocity  components  in  directions  respectively 


- - Shear  stress 
P - Field  density 
y - Stream  function 
® - Electrical  conductivity 
K -consistency  coefficient 
' - Shear  rate 

Bu 

u - Uniform  strength  of  magnetic  Jield 
^ - Power-Iaw  index 

n 

1 - Similarity  variab1e 
^ - Velocity  parameter 
M . Magnetic  parameter 
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1.  INTRODUCTION 


Researchers  in  mechanical  discipline  have  curiosity  and  interest  to  find  the  solution  of  nonlinear  differential 
equations  for  governing  mathematical  model.  Applications  of  boundary  layer  flow  are  in  textile  and  paper  industry, 
manufacture  of  sheets,  crystalline  materials,  spinning  of  fibers.  [l]Wu  (1973)  has  studied  the  effects  of  suction  or  injection 
on  MHD  boundary  layer  flow.  [2]  Takhar  et.al.  (1987)  studied  a MHD  asymmetric  flow  over  a semi-infinite  moving 
surface  [3]  Mahapatra  and  Gupta  (2001)  studied  a steady  two-  dimensional  stagnation-point  flow  by  a uniform  transverse 
magnetic  field.  [4]  Jean-David  Hoernel  (2008)  has  been  investigated  the  similarity  solutions  for  the  steady  laminar 
incompressible  boundary  layer  governing  MHD  llow.  [5]  Rajput  (2014)  study  of  MHD  boundary  layer  flow  of  non- 
newtonian  power  law  fluid.  Pakdemirli  [6]  derived  the  boundary  layer  equations  of  power-fluids.  Cubic  Spline  functions 
were  used  to  solve  second  order  Foppl-Hencky  equation  by  [7]  Pandya  et  al.  (2010). 

The  object  of  the  present  paper  is  to  study  the  effect  of  magnetic  parameter  on  MHD  flow  for  two  cases,  when 
plate  is  at  rest  and  lluid  and  plate  moves  in  same  direction  with  same  velocity.  The  governing  equation  is  nonlinear 
differential  equations,  which  is  solved  by  using  spline  collocation  method  due  to  blue  [8].  In  this  way,  the  paper  has  been 
organized  as  follows.  In  section  2,  problem  formulation  is  given,  in  section  3,  method  description  and  approximate  solution 
by  using  collocation  method  and  results  and  discussion  are  presented  in  section  4.  Section  5 includes  the  conclusion. 

2.  GOYERNING  EQUATIONS 


Here  we  studied  the  two-dimensional  laminar  boundary  layer  flow  of  a viscous,  incompressible  and  electro 
conducting  power  law  fluid  past  a continuously  moving  surface  passing  through  with  constant  w in  the  same  direction  to 


the  free  stream  velocity  00 . The  - axis  extends  parallel  to  the  plate  and  -y  - axis  perpendicular  to  the  A - axis.  A 

magnetic  field  of  uniform  strength  js  appjjec|  jn  the  positive  - direction,  which  produce  the  magnetic  field  in  the  x - 
direction.  The  boundary  layer  equations  governing  the  flow  in  a power-law  fluid  are 


y. 


X 


dii  dv  _ 
r)x  c)y 


(2.1) 


c)u  dv  1 du  (j  Br. 

u — + v — = -u 

dx  dy  p dy  p 


(2.2) 


Where  ^ are  the  velocity  components  along  and  coordinates,  xv  is  the  shear  stress  and  P is  the  lluid 

density. 

with  the  boundary  conditions: 


y = 0:  U = UW,  v = 0 
y = oo : u = Um 

We  apply  power-law  relation  between  the  shear  stress  and  the  shear  rate  by 


(2.3) 
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-JjL-c  = 

dx 


V- : _T  = K 

v — _ xy 


du 


dy 


n—1 


du 

dy 


r 


Where 


dlt 


dy 


n- 1 


r- 


K 


denotes  the  kinematic  viscosity,  K is  the  consistency  coefficient  P anc)  n js  the  power- 


law  index,  for  n<  \ pSeucj0  plastic,  for  n ^ the  fluid  is  Newtonian,  n ^ ^ for  dilatants  fluid.  The  equation  (8)  becomes 


du  3v  d 

u h V — = — 

c)x  dy  dy 


r 


du 


dy 


n- 1 _ A 


du 

dy 


<jB- 


(2.4) 


W(x,  y) 


by/  dy/ 
u = — y = — — 

dx 


Introducing  the  stream  function  + ’ ■■  ’ such  that  - , (,x  , which  satisfy  the  continuity  equation 

(2.1).  Govind  R.  Rajput  et  al  [5]  converted  partial  differential  equation  into  nonlinear  ordinary  differential  equations  using 
Group  theoretic  method.  They  considered  the  following  transformation: 


y/{x,  y)  = axaf(7]),  ij  = b\ 


(2.5) 


Where  a^bW/  anc|  P are  real  numbers,  ^ is  similarity  variable,  f ■ n ^ is  the  transformed  dimensionless  stream 


function. 


Applying  this  similarity  variable  ^ they  derive 

¥x  = axa  ] [a  f - Tjfif  j 
y/y  = abf'xa~p 
Wyy  =ab2xa~2fS  f 

¥yx  = abxa-p~\af'  -Pf'  (2  6) 

Using  the  equation  (2.5)  along  with  (2.6)  into  (2.4)  they  get  transformed  into  nonlinear  ordinary  differential  of  the 

form 


(l/T  / j ~Mf +di#" =0  (l/'r  /') -Mf+yrr =0 

With  the  transformed  boundary  conditions: 

/ (0)  = 0,  /(0)=G,/(oo)=l 

n =1  gives  Newtonian  fluids,  then  the  equation  (2.7)  becomes 
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With  boundary  conditions 
/<0)  = 0,  / (0)  =e,  / (oo)  = l 

S=T 

Where  “ is  the  velocity  parameter  and 


(2.9) 


(2.10) 


M = ^ 


pu„ 


is  the  magnetic  parameter.  Here  note  that  when 
^ — 0 plate  is  stationary  and  e = 1 plate  and  fluid  moves  same  direction  and  same  velocity. 

3.  QUARTIC  SPLINE  BLUE  METHOD 

For  three  points  boundary  value  problems  are 

s.  (x)  [ X, , X 1 

Let  ' be  quartic  spline  runction  in L J . Conditions  for  natural  splines  are 

s.  (x)  [ x_, , X 1 

• ' Almost  quartic  tn  each  subinterval L J . 

. •s/(-*i)  = yi  , for  i = 0,1,2, ,n. 

.s-;(x;)  , S,.(x,.),  S,.  (x; ),  s (x;.)  are  continuous  in  [x0xn]. 


s,'(x o)  = s:\xn)  = 0. 

[ x_, , x 1 

Here  spline  third  derivative  must  be  linear  in L J . So, 

V(*i)  = H (*,  - x)yf +(x-x!-i)y'] 


Where  = X'  _ ^-1  and  s'Uii  = yi 


Integrate  (3.1)  twice  with  respect  to  X 


si(x)  = 

h. 


+ c,  (x;.  - x)  + r/,  (x  - x;_j ). 


(3.1) 


Where  use 


A'U,-i)  = .Vi  and si(xi)  = yi  i 


c.  . d. 


1 

c.  = — 


h. 


h; 


'i  V 


Jm  - — 1 

6 


2 and 


J,=- 

/r 


in  (3.1),  we  get  constants  ' and  ' 


2 A 


li  V 


• \ 

yt  -—yt 

o 
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So 


*,•'(*)  = y 

k; 


r y t— i ^ z 


i 

+— 

/t 


/i.2 


y,--i  -^>,-1 

v 6 y 


(*,-  -+)  + 


l 


/i. 


2 A 


■ h., 

y,  ~—yi 


'i  V 


y 


Integrate  (3.2),  once  with  respect  to  x, 

1 


5,U)  = 


h 


y,-t  + „ „ y,- 


24 


24 


h 


hr  ,„l(x-x)2 


t V 


y,-,  -— y,-i 

o 


+ 


1 

+— 


h 


h: 


7 V 


yi 

0 


(x- x,_|)2 


Take 


S,  (X_, ) = V.  , e, 

1 11  1 1 , we  get  constants  1 


h3  . /7  . 

where  ° 4 


e 

Substitute  ! in  (3.3),  we  get 


•>(+)  = — 
h 


ixi~x ) • {x~xi- 1)  ■ 

>i-i  + „ „ >,■ 


24 


24 


/7. 


V 


O 


(>,->)  2 


+ 


1 

+ — 


/7,- 


i V 


>,  “ — >Z 
o 


(jc-x_,)2  h3  . /7 

+ >/-|-T>,-1  +->,-,• 


Here 


,v(  (x;  ) = si+1  (x;+)  and  for  equal  intervals  we  have, 


, . h2  ,,,  

>/+l  “2>/  +>,_!  =—(>/+i  +4 >,-  +>,_!  ) 

o 


And  for  'S'^X'  ) = ‘s'<+,  (>,+ ) and  for  equal  intervals  we  have, 
h ( , . \ h3 

>,-,  ->,  =--(>,  + >/-, ) + — (3>,-,  - >/  ) 

3.1  Solution  of  the  Problem  by  using  Spline  Function  Due  to  Blue 

We  solve  the  nonlinear  ordinary  differential  equation 


(3.2) 


(3.3) 


(3.4) 


(3.5) 


(3.6) 
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f-W+\ff'=  0 

With  boundary  conditions 

/ (0)  = 0,  /'(0)=e,/  (l)  = l 


(3.1.1) 


(3.1.2) 


Solve  this  for  two  cases  when  ^ ® plate  is  stationary  and  £ — 1 plate  and  fluid  moves  same  direction  and  same 

velocity. 

Case  (i)e  — ® plate  is  stationary 

/ (V  ) = 2 

To  obtain  the  spline  solution,  we  begin  with  an  assumed  function  4 , which  satisfy  given  boundary 

conditions  (3.1.2).  To  find  the  solution  of  equation  (3.1.1)  along  with  boundary  conditions  (3.1.2),  we  use  2 , 

V-  / = 1214 

equations  (3.1.1)  and  (3.5).  Put  h=0.1,  we  get  different  values  of  • ' for  ’ ’ ’ ’ 

• |234 

To  find  the  final  solution  we  use  (3.6)  for  different  values  ol'  1 ~ ’ ’ ’ 4 respectively 
Thus,  we  obtain  the  system  of  equations  as  follows 

/o  - /i  = + /oJ^L-t,"+3To"J 

h . . h 3 

/l-/2  =--[/2  +/J  + — L-/2  +3/l  J 

h . l f . ...  . » . 

/2-/3  = -^[/3  +/2J  + — [-/3  +3/2  ] 

h . //  

/3 -/4  = --[/4  +y3]  +—i-y4  +3y3  ] 

2 24  (3.1.3) 


Substitute  and  tor  ^ L).l 
Matlab,  we  get  solution  graph  as  follows: 


in  (3.1.3).  We  get  four  unknown  and  four  equations.  Solving  that  system  using 
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Plat  stationary 


Figure  1:  /3//  ■ versus  ^ for  different  Values  of  M 

Plat  stationary 


f (n)  77 

Figure  2:  ■’  u versus  ' for  different  Values  of  M 
Case  (ii)G  = 1 plate  and  Huid  moves  same  direction  and  same  velocity 

To  obtain  the  spline  solution,  we  begin  with  an  assumed  function  / ~ which  satisfy  given  boundary 

conditions  (3.1.2).  To  find  the  solution  of  equation  (3.1.1)  along  with  boundary  conditions  (3.1.2).  First  we  use  f^^  ~ ^ , 

equation  (3.1.1)  and  (3.5)  and  0.1  to  get  different  values  of  for  ; — 1>2,3,4. 


To  find  the  final  solution  we  use  (3.6)  for  different  values  of  1 3, 4 respectivcly  and  get  the  system  of 


equations 


h , , h 3 

To-Ti  = --[?!  +y0}+-\-yi  +3t0  ] 

k r ■ >,  /t3  r - 0 

+yJ+— [-y2  +3^] 

h , , /r’  

yi-y3  = --ly3  +t2]+— [-^  +3t2] 
y3  - = ~[y4  + y3] + -y*  + 3 y3'] 

z (3.1.4) 


Substitute 


y'i 


and 


yt 


for /1  = 0.1 


in  (3.1.4).  We  get  four  unknown  and  four  equations,  solving  which  using 
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Matlab  we  get  solution  graph  as  follows: 


Fluid  and  Plate  moves  same  direction 


Figure  3: 


/(7) 


versus 


7 


for  different  values  of  M 


Fluid  and  Plate  moves  same  direction 


Figure  4:  ' ^ versus  ^ for  different  values  of  M 

Here  we  have  shown  the  comparison  of  velocity  and  displacement  when  plate  is  stationary  and  plate  and  fluid 
moves  in  same  direction  with  same  velocity. 
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Figure  5:  Comparison  of  $ Yersus^  , when  M=1 
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1 


Figure  6:  Comparison  of  ' ^ Yersus  ^ , when  M=1 

4.  RESULTS  AND  DISCUSSIONS 

The  result  shows  that,  there  is  a significant  impact  of  magnetic  parameter  on  the  displacement  profile  of  the  flow 
(Figure  1,  3).  In  both  the  cases  we  can  see  that,  displacement  of  fluid  increase  exponentially  but  as  magnetic  parameter 
increases,  displacement  of  fluid  decreases.  Figure  2 demonstrate  that  velocity  of  MHD  power  law  Newtonian  fluid  is 
transversely  proportional  to  magnetic  parameter.  When  fluid  and  plate  moves  in  same  direction  with  same  velocity,  there  is 
no  change  in  velocity  but  increment  in  magnetic  parameter  shows  remarkable  change.  Velocity  profile  undergoes  decay  as 
increase  in  magnetic  parameter  Figure  4.  From  figure  5,  when  plate  is  stationary,  displacement  of  lluid  is  lesser  than 
displacement  of  fluid  when  plate  and  fluid  moves  same  direction.  Comparisons  of  velocity  profile  in  both  the  cases  are 
given  in  figure  6. 

5.  CONCLUSIONS 

In  power-law  MHD  flow,  there  is  remarkable  effect  of  magnetic  parameter  on  flow,  whether  the  plate  is  at  rest  or 
Huid  and  plate  moves.  Magnetic  parameter  is  transversely  proportional  to  the  velocity  of  fluid. 

We  find  the  generalization  of  blue  method  for  third  order  problem  and  solved  the  problem  using  spline  collocation 
method  due  to  Blue.  The  beauty  of  this  method  is,  there  is  no  need  to  convert  nonlinear  problem  into  linear  one,  and  we 
can  get  the  solution  directly  for  the  nonlinear  form.  Thus,  we  can  get  the  results  for  practical  approach,  without  doing  the 
experimental  work  and  save  the  time  and  experimental  cost. 
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